Search for: All records

Non-adiabatic (NA) molecular dynamics (MD) is a powerful approach for studying far-from-equilibrium quantum dynamics in photophysical and photochemical systems. Most NA-MD methods are developed and tested with few-state models, and their validity with complex systems involving many states is not well studied. By modeling intraband equilibration and interband recombination of charge carriers in MoS2, we investigate the convergence of three popular NA-MD algorithms, fewest switches surface hopping (FSSH), global flux surface hopping (GFSH), and decoherence induced surface hopping (DISH) with the number of states. Only the standard DISH algorithm converges with the number of states and produces Boltzmann equilibrium. Unitary propagation of the wave function in FSSH and GFSH violates the Boltzmann distribution, leads to internal inconsistency between time-dependent Schrödinger equation state populations and trajectory counts, and produces non-convergent results. Introducing decoherence in FSSH and GFSH by collapsing the wave function fixes these problems. The simplified version of DISH that omits projecting out the occupied state and is applicable to few-state systems also causes problems when the number of states is increased. We discuss the algorithmic application of wave function collapse and Boltzmann detailed balance and provide detailed FSSH, GFSH, and DISH flow charts. The use of convergent NA-MD methods is highly important for modeling complicated quantum processes involving multiple states. Our findings provide the basis for investigating quantum dynamics in realistic complex systems. 

Nonadiabatic molecular dynamics (NA-MD) is a powerful tool to model far-from-equilibrium processes, such as photochemical reactions and charge transport. NA-MD application to condensed phase has drawn tremendous attention recently for development of next-generation energy and optoelectronic materials. Studies of condensed matter allow one to employ efficient computational tools, such as density functional theory (DFT) and classical path approximation (CPA). Still, system size and simulation timescale are strongly limited by costly ab initio calculations of electronic energies, forces, and NA couplings. We resolve the limitations by developing a fully machine learning (ML) approach in which all the above properties are obtained using neural networks based on local descriptors. The ML models correlate the target properties for NA-MD, implemented with DFT and CPA, directly to the system structure. Trained on small systems, the neural networks are applied to large systems and long timescales, extending NA-MD capabilities by orders of magnitude. We demonstrate the approach with dependence of charge trapping and recombination on defect concentration in MoS₂. Defects provide the main mechanism of charge losses, resulting in performance degradation. Charge trapping slows with decreasing defect concentration; however, recombination exhibits complex dependence, conditional on whether it occurs between free or trapped charges, and relative concentrations of carriers and defects. Delocalized shallow traps can become localized with increasing temperature, changing trapping and recombination behavior. Completely based on ML, the approach bridges the gap between theoretical models and realistic experimental conditions and enables NA-MD on thousand-atom systems and many nanoseconds. 

Grain boundaries (GBs) in perovskite solar cells and optoelectronic devices are widely regarded as detrimental defects that accelerate charge and energy losses through nonradiative carrier trapping and recombination, but the mechanism is still under debate owing to the diversity of GB configurations and behaviors. We combine ab initio electronic structure and machine learning force field to investigate evolution of the geometric and electronic structure of a CsPbBr 3 GB on a nanosecond timescale, which is comparable with the carrier recombination time. We demonstrate that the GB slides spontaneously within a few picoseconds increasing the band gap. Subsequent structural oscillations dynamically produce midgap trap states through Pb–Pb interactions across the GB. After several hundred picoseconds, structural distortions start to occur, increasing the occurrence of deep midgap states. We identify a distinct correlation of the average Pb–Pb distance and fluctuations in the ion coordination numbers with the appearance of the midgap states. Suppressing GB distortions through annealing and breaking up Pb–Pb dimers by passivation can efficiently alleviate the detrimental effects of GBs in perovskites. The study provides new insights into passivation of the detrimental GB defects, and demonstrates that structural and charge carrier dynamics in perovskites are intimately coupled. 

Photoinduced nonequilibrium processes in nanoscale materials play key roles in photovoltaic and photocatalytic applications. This review summarizes recent theoretical investigations of excited state dynamics in metal halide perovskites (MHPs), carried out using a state-of-the-art methodology combining nonadiabatic molecular dynamics with real-time time-dependent density functional theory. The simulations allow one to study evolution of charge carriers at the ab initio level and in the time-domain, in direct connection with time-resolved spectroscopy experiments. Eliminating the need for the common approximations, such as harmonic phonons, a choice of the reaction coordinate, weak electron–phonon coupling, a particular kinetic mechanism, and perturbative calculation of rate constants, we model full-dimensional quantum dynamics of electrons coupled to semiclassical vibrations. We study realistic aspects of material composition and structure and their influence on various nonequilibrium processes, including nonradiative trapping and relaxation of charge carriers, hot carrier cooling and luminescence, Auger-type charge–charge scattering, multiple excitons generation and recombination, charge and energy transfer between donor and acceptor materials, and charge recombination inside individual materials and across donor/acceptor interfaces. These phenomena are illustrated with representative materials and interfaces. Focus is placed on response to external perturbations, formation of point defects and their passivation, mixed stoichiometries, dopants, grain boundaries, and interfaces of MHPs with charge transport layers, and quantum confinement. In addition to bulk materials, perovskite quantum dots and 2D perovskites with different layer and spacer cation structures, edge passivation, and dielectric screening are discussed. The atomistic insights into excited state dynamics under realistic conditions provide the fundamental understanding needed for design of advanced solar energy and optoelectronic devices.